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ABSTRACT 

We study the spectral energy distribution of gamma rays and neutrinos in the process- 
ing microquasar SS433 as a result of pp interactions within its dark jets. Gamma-ray 
absorption due to interactions with matter of the extended disk and of the star is 
found to be important, as well as absorption caused by the UV and mid-IR radia- 
tion from the equatorial envelopment. We analyze the range of precessional phases for 
which this attenuation is at a minimum and the chances for detection of a gamma-ray 
signal are enhanced. The power of relativistic protons in the jets, a free parameter 
of the model, is constrained by HEGRA data. This imposes limits on the gamma-ray 
fluxes to be detected with instruments such as GLAST, VERITAS and MAGIC II. 
A future detection of high energy neutrinos with cubic kilometer telescopes such as 
IceCube would also yield important information about acceleration mechanisms that 
may take place in the dark jets. Overall, the determination of the ratio of gamma-ray 
to neutrino flux will result in a key observational tool to clarify the physics of heavy 
jets. 
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1 INTRODUCTION 

The famous and enigmatic microquasar SS433 has been mat- 
ter of investigation for more than two decades. Consisting 
of a donor star feeding mass to a black hole, it presents two 
oppositely directed, precessing jets with hadronic contenlQ. 
We refer to the rel ativistic collimated outflows as 'dark' jets 
l|Gallo et all 120051 ) since the very hi gh kinetic luminosity 
~ 10 39 erg s _1 (|Dubner et al.|[l998h appears as the domi- 
nant power output of the ejected material, having imprinted 
a deformation on the supernova remnant W50. 

Most of the radiative output of the system is observed 
in the UV and optical bands, whereas the X-ray emission 
detected is ~ 1000 lower than the kinetic energy of the 
jets, pro bably due to a screening effect with an equa torial 
outflow (|Revnivtsev et al.ll2006l ; Marshall et al.ll2006h . The 
gamma-ray emission above 0.8 TeV has be en constrained by 
HEGRA to be $ 7 < 8.93 x lO^cm"^ -1 i|Aharonian et all 
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1 Iron lines with a shift corresponding to a velocity of v ~ 0.26c 
have been detected, for instance, by Migliari et al. (2002). 



2005) whereas the neutrino flux upper limit accor ding to 
AMA NDA-II data is < 0.21 x lO^cm^s -1 (|Halzenl 
l2006h . 

In previous hadronic models for high energy emission 
from microquasars, relativistic protons in the jets inter- 
act with tar get protons from the stellar wind of the com- 
panion star llRomero et alj 120031 : IChristiansen et all 120061 : 
lOrellana et al.ll2007l h Since in the case of SS433 there is no 
evidence of such a strong stellar wind, in this work we inves- 
tigate the possible production of gamma rays and neutrinos 
resulting from pp interactions between relativistic and cold 
protons within the jets themselves. 



2 PRELIMINARIES 

The binary SS433, distant 5.5 kpc from the Earth, dis- 
plays two mildly relativistic jets (vb ~ 0.26c) that are op- 
positely directed and precess in cones of half opening an- 
gles of 6 « 21°. The line of sight makes an angle i = 78° 
with the normal to the orbital plane and a time-dependent 
angle ij(t) with the approaching jet (see Fig. [TJ. Assum- 
ing that ip(t) is the precessional phase of the approach- 
ing jet, we shall follow the convention that when ip = 
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approaching jet 



extended disk 




where the kinetic energy density is 



receding jet 



Figure 1. Schematic view of the SS433. The approaching jet is 
most of the time closest to our line of sight and the receding one 
is oppositely directed. 



the mentioned jet points closer to the Earth. Then, when 
ip = 0.5, it has performed half of the precession cycle and 
it makes its largest angle with the line of sight. The mass 
loss rate in the jets is rhj = 5 x 10~ Moyr -1 , the pe- 
rio d of precession is 162 d and the orbital period is 13.1 
d l|Fabrikal 120041 ') . The donor star and the compact ob- 
ject are thought to be embedded in a thick expanding disk 
which is fed by a wind f rom the supercritica l accretion disk 
around the black hole l|Zwitter et al.| [T991). This equato- 
rial env e lope i s perpendicular to the jets and according to 
iFabrikal l|2004h we assume that it has a half opening angle 
a w « 30°, a mass loss rate M w ~ W~ 4 Mqyv^ 1 and a ter- 
minal velocity v v ~ 1500 km s _ . Also, this extended disk 
has been rec ognized as the origin of both the UV and mid- 
IR emission ijGies et al.ll2002al : lFuchs et al.ll2005h which can 
cause significant abs orption of gamma-rays as discussed in 
iRevnoso et ail (|2008T ). 

The spectral identification of the companion star has 
been difficult due to the presence of the extended disk, since 
the star is often partially or totally obscured by it. After 
convenient observations at specific configurations of preces- 
sional and orbital phases it has became quite clear that the 
star i s an A-supergiant dHillwig et al.l |2004| ; iBarnes et al.l 
l200rj ; ICherepashchuk et al.ll2005l ). We assume the masses of 
the components as derived from INTEGRAL observations 
jCherepashchuk et alJliooBI ). M hh = 9M and Af* = 3OM 
for the black hole and star respectively. This corresponds to 
an orbital separation a ~ 79 Rq for a zero-eccentricity orbit 
as it is the case for SS433. Since the star is believed to fill it s 
Roche lobe, the implied radius according to lEggletonl (|l983l ) 
is R L ~ 38R Q . 



2.1 Outline of the jet model 

We assume that a magneto-hydrodynamic mechanism for 
jet ejection operates in SS433, that is, ejection is realized 
through the conversion of magnetic energy into matter ki- 
netic energy. The magnetic energy density is supposed to 
be in equipartition with the kinetic energy density of the 
ejected particles, so that the corresponding magnetic field 
along the jet is given by 



m p v b -n-R? (zj) 



(2) 



Here E^ is the classical kinetic energy of a jet proton with 
velocity v b and Rj(z) is the jet radius at the height Zj along 
the jet axis. 

The j ets are modeled as co nes with a half opening angle 
£j ~ 0.6° (Ma rshall et alj [20061 ). Assuming an initial jet ra- 
dius Ro — Rj(zo) ~ 5J?sch, where iisch = 2GMbh/c 2 , we find 
the injection point as zn = i?o/tan£j ~ 1.3 x 10 9 cm. Since 
the jets are heavy as compared to other similar objects, it is 
reasonable to admit that they are cold matter dominated. In 
this case, we assume that a small fraction of relativistic or 
hot particles are confined by the cold plasma. According to 
iBosch-Ramon et al.l (|2006l ) the pressure of cold particles is 
greater than that of the relativistic ones if the ratio of cold 
to hot particles is less than 1/1000, and this condition will 
be greatly satisfied provided that the luminosity carried by 
relativistic particles is required to be smaller than the total 
kinetic luminosity of the jet. 

Particle acceleration is supposed to take place via dif- 
fusive acceleration by internal shocks converting bulk ki- 
netic energy into random kinetic energy. According to the 
standard model for non-relativistic shock acceleration (e.g. 
Blandford & Eichler 1987 and references therein) we expect 
that the relativistic proton spectrum is given by a power-law, 
Np(E' p ) — KoE'~ a at zj = Zq, where the spectral index is 
the standard value for first order diffusive shock accelera- 
tion, a = 2. The flux of these protons hence evolves with Zj 
as 

Jp( E 'p) = -^-( z o/zj) 2 E' p a 

in the jet frame, which transform ed to the observer frame 
l|Purmohammad fc Sa mimi 200l]) reads 



Jp{t> E p , Zj 



cK ( z 

47T 



r- Q+1 (E P -p b ^Ef 



rripC 4 cos ij) 



sin 2 ij + T 2 cosij(t) — 



l3 b E p 



Z() 



J P {Ep,t), (3) 



where ij (i) is the angle between the jet axis and the line of 
sight, /3b = 0.26, and T = [l — /3 2 ] ^ is the jet Lorentz 
factor. The normalization constant Ko is obtained by speci- 
fying the fraction of power carried by the relativistic protons, 

9rel, 

,/(max) 

(4) 



/•Bp 

ttRI / Jp{E'p)E'pdE'p = q re \L k 



IE 

so that 



A', 



o 



cRl In 



E'(, min > 
/(min) 



(•») 



1 GeV and the maximum proton 



(1) 



where we take E 
energy E' p ma ^ will be determined in the next section. We 
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shall adopt, for the illustrative predictions of neutrino and 
gamma-ray fluxes, a tentative value q Te i = 10 -4 , but a full 
discussion of the possible range for this parameter will be 
presented in Sect. [5] 



3 HADRONIC PROCESSES IN THE JETS 

Relativistic protons in the jets are subject to different mech- 
anisms that can make them lose energy. In this section we 
analyze the energy range where pp collisions are the dom- 
inant cooling process that will produce the corresponding 
gamma rays and neutrinos in SS433. 



3.1 Acceleration 

The accelerat ion rate of protons up to an energy E p can be 
estimated as (|Begelman et al.lll990l ): 



t: 



ceB 
E„ ' 



(6) 



where r\ ~ f3 h is the acceleration efficiency. Clearly, as long 
as the latter rate is greater than the total loss rate for a 
given energy, protons will be effectively accelerated up to 
that energy. 

As mentioned in Sect. 2.1, we assume that the protons 
are accelerated at shocks produced by collisions of plasma 
outflows with different bulk velocities. In the frame of the 
shock, the conservation equations imply that the upstream 
velocity is significantly higher than the downstream veloc- 
ity, i.e., v u /vd = (Theat + l)/(7hcat — 1) = £, with 7 hca t the 
ratio of specific heats and £ the compression factor. In this 
way, the two regions may be regarded as two converging 
flows. The Fermi first order acceleration mechanism then 
operates to produce a power law particle spectrum, which is 
essentially independent of the microphysics involved. Strong 
shocks (£ ~ 4) can be non-relativistic as it is the case, for 
instance, in supernova remnants and colliding wind massive 
binaries. In the case of sub-relativistic jets, strong shocks are 
expected as suggested by the non-thermal synchrotron ra- 
dio spectra observed from the jets of microquasars (Fender 
2004). Non- linear effects like shock modification by the pres- 
sure of the relativistic particles or magnetic field effects can 
result in a variety of spectral indexes. The reader is referred 
to the recent paper by Rieger et al. (2006) on Fermi accelera- 
tion in astrophysical jets, which includes a section on mildly 
relativistic microquasar outflows. 

The acceleration of protons proceeds only for protons 
with a threshold energy that allows the diffusive accelera- 
tion process to take place (Rieger et al. 2006). Then, only 
the supra-thermal tail of the Maxwellian distribution of cold 
protons will be affected by the process (Bosch- Ramon et al. 
2006). This has the result that just a small fraction of the 
total power carried by the jet is converted to relativistic 
particles. 

3.2 Cooling rates and maximum particle energy 

The density of cold protons at a distance Zj from the black 
hole in each jet is 



n p ( Zj ) 



■K[Rj(zj)] 2 m p v h ' 



(7) 



These cold protons serve as targets for the relativistic ones, 
so that the rate of pp collisions in the jet is given by 



(8) 



where the inelasticity coefficient is taken to be K pp ~ 1/2 
since on average, the leading proton losses half of its total 
energy per collision. 

The cross se ction for inelastic p p interactions can be 
approximated by (|Kelner et al.|[20060 



j pp \e p ) 



(34.3 + 1.88L + 0.25L 2 



Eth 

E V 



, -in-27 2 

x 10 cm , 



(9) 
(10) 



where L = ln(£ p /1000 GeV) and E th = 1.22 GeV. 

Cooling by jry interactions can take place mainly via 
photomeson pr oduction (jp — > pn % ) and pair production 
(jP —* pee + ) (Be gelman et al.l 119901). The corresponding 
cooling rate can be obtained from ( Atovan fc Dermer]l2003 ) 



dE 



cn ph (E) 



2l p E 



(11) 



where E{. h ~ 150 MeV, y p is the Lorentz factor of the proton, 
n p h{E) represents the density of target photons, o vl will 
be the inelastic cross section appropriate for photopion and 
photopair creation, and K pl is the corresponding inelasticity 
coefficient. 

Photopion production will occur when protons 
collide with X-ray phot ons, for which, based on 
Chcrcpashchu k~et al.l l|2005h . we adopt a Bremsstrahlung 
X-ray distribution for 2 keV < E < 100 keV, 



nx(-B) = Lx- 



-E/(kT e ) 

4tvz?E' 2 



(erg 'cm 



(12) 



^ 36 erg s _1 . These X- 



where kT e « 30 keV and L x = 
ray photons are considered to be originated in a corona 
surrounding the inner ac cretion disk, as suggested in 
ICherepashchuk et all (|2005h . 

The c ross section for photopio n production is approxi- 
mated by (jAtovan fc Dermer|[2003l l 

o$ = e(£ r -200 MeV)6(500 MeV-£ r ) 3.4 x 10~ 28 cm 2 
+ Q{E r - 500 MeV) 1.2 x 10" 28 cm 2 , (13) 

and the inelasticity coefficient for photopion production is 



K 



<d(E r - 200 MeV)e(500 MeV - E r ) 0.2 

+ Q(E r - 500 MeV) 0.6. (14) 



The contribution of the e~e pair creation process to 
the total p7 cooling rate is calculated also using Eq. (JTTJ , but 
the soft photon density in this case includes also the contri- 
bution associated with the UV emission from the extended 
disk, 



n ph (E) = n vv (E,n) 



+ n x (E). 



Based on the discussion in iGies et al.l l|2002ah , the UV 
photons with wavelengths in the range (1000 A, 10000 A) 
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correspond to a blackbody distribution with Tuv = 21000 
K. Hence, we take the corresponding radiation density as 



nvv(E,Q,) = 2E 2 (hc) 



-3,E/kT uv 



(15) 



For this process we consid er the expressions for cross 
section and inelasticity given in Begel man et all ljl990h : 

<r$ = 5.8 x 10~ 28 cm 2 [3.111/ - 8.07+ 

(2m e c 2 /E) 2 (2.7 X' - L' 2 + 0.67L' 3 + 0.55) - 
(2m e c 2 /E) i (0.19L' + 0.13) - (2m e c 2 /E) 6 (0.01L')] 



and 



K w = 4- 

m v E r 



-8.78 + 5.51 V - 1.61 L' 2 + 0.69 L' : 



3.11 L' - 8.07 



(16) 



(17) 



with L' = In { 

The accelerated protons can also lose energy in the form 
of synchrotron radiation at a rate 



a T B 2 
m e c 8ir 



7p, 



(18) 



and via Inverse Compton scatterings with X-ray and UV 
photons at a rate 



m e c 



7p, 



(19) 



where 



e p h 



'/Bp 



n ph (_B) E dE 



E min 



is the corresponding d ensity of energy in target soft photons 
jBegelman et alj|l990h . 

It is also expected that the accelerated protons will suf- 
fer adiabatic losses because of the expansion undergone by 
the jets. The corresponding adiabatic cooling rate can be 
written as (see e.g. Bosch Ramon et al. 2006) 



t: 



2«b 
3 Zj 



(20) 



We show the obtained results for the acceleration and 
cooling rates at the base of each jet (zj — Zq) in Fig. [5] 

As it can be seen from this plot, the expected cooling 
rates for p^y and Inverse Compton interactions are found 
to be negligible. We infer that the protons which can be 
effectively accelerated up to energies below 

(max) 



E 



3.4 x 10 D GeV, 



(21) 



will cool efficiently mainly via pp collisions at Zq. For higher 
energies, synchrotron losses become dominant. On the other 
hand, as z, increases, adiabatic losses also grow and the max- 
imum proton energy at which pp collisions dominate can 
change along the jets. 

The absolute maximal energy for protons at a given Zj, 
X p abs) (see Fig. |3J) is obtained from t^ ccl = + t^ lah + 
t~ l 4- r _1 4- r _1 

''sync ' L IC ^ VI ' 

The size constraint, implying that the proton gyro- 
radio has to be smaller than the radius of the jet, i.e. 



E p < Ep 1 



(size) 



eRjB ~ 3 x 10 GeV, does not happen to 
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limit the energy of the protons at the bases of the jets. Note 



10 1 10 2 10 3 10 4 10 5 10 6 10 7 10 s 
E p [GeV] 

Figure 2. Proton accelerating and cooling rates at the base of 
the jets. 



also that we have the same value of _Ep Slzc ' for larger values 
of Zj along the jets because B oc RT 1 . Nevertheless, the size 
constraint will limit the energy of the accelerated protons at 
distances Zj > 3 x 10 12 cm from the black hole (see Fig. [3). 

In the (zj,E p ) region where pp collisions dominate the 
cooling mechanism, the condition t~p > t'y 1 ^ +<^d la b+iic + 
tp} must hold. This region is indicated in the shaded zone of 
Fig. [3] It can be seen that the maximum energies for efficient 
cooling through pp interactions are Ep ~ 3 x 10 6 GeV for 
Zj < 0i, where Zj ~ lO 12 cm0 Therefore, the jet will become 
essentially cold and observationally dark, unless some re- 
acceleration mechanism could operate (e.g. mediated by re- 
collimation shocks or terminal shocks). 

It is also worth noting that decreasing the parameter 
q Te i does not change the maximum proton energy given by 
()21l) since in that case the target proton density is essentially 
unchanged and so is the pp cooling rate. This allows us to 
take g rc i as a free parameter that can be easily factored out 
in all of our predictions for gamma-ray and neutrino signals. 



3.3 High-energy gamma rays from pp interactions 

The collision of a certain number of high-energy protons 
with cold protons in the jets will cause the production of 
second ary gamma rays and neutrinos. Following the treat- 
ment of lKelner et~all (|2006l ). which is based on SIBYLL sim- 
ulations of pp interactions including perturbative QCD ef- 
fects, the spectrum of produced gamma-rays with energy 



2 Since the necessary maximal energy in equation 10 has to be 
expressed in the frame comoving with the jet, we take E'( > max ) ~ 

- m^c 4 ), where _B^, max) ps 3.4 x 10 6 GeV 



/W4 max) " 

at the base of the jet. 
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Figure 3. The region in the (zj,E p ) plane where pp interac- 
tions are the dominant cooling mechanism is indicated in gray. 
The absolute maximum proton energy (solid line) and the maxi- 
mum proton energy from the size constraint (dotted line) are also 
shown. 




Figure 4. Spectral intensity of gamma-rays produced in the ap- 
proaching jet as a function of the precessional phase and the 
gamma-ray energy. 



where E m in = E-, + ^f^, E max = K 1T (E { p nz ' x ' ) -m p c 2 ), and 



E-y — xE p for a primary proton with energy E p reads 



F^(x,E p ) = B. 



\nx 



,/3 7 



1 

In x 



1 + fc 7 aA(l - x p i) 
4/3 7 ^ 4fc 7 /3 7 x^(l - 2x f3 ~>) 



■x p i 1 + kyX p i (1 - xP-t) 



where 
By = 

Pj = 



1.3 + 0.14 L + 0.011 L 
1 

1.79 + 0.11 L + 0.008 L 2 
1 

0.801 +0.049 L + 0.014 L 2 ' 



(22) 

(23) 
(24) 

(25) 



with L = In (E p /1 TeV). 

For _E 7 > 100 GeV, we shall consider the gamma-ray 
emissivity at a height zj along the jets as 



dN~,(t,E.y, Zj) 
d~K 



E*y 



J' 



r [x,^- \ <!.<■ (2()i 



z\ I dE~ 



(27) 



in units GeV _1 s _1 . The integration limits Xmin and s max 
are chosen in order to cover the proton energy range where 
pp collisions dominate at each Zj, as shown in Fig. [3] On 
the other ha nd, for _E 7 < 100 G eV, we shall consider, as 
suggested in iKelner et alj l|2006l ), the emissivity obtained 
using the (5-functional approximation 



qn(t,E n ,Zi 



n incl/ 2 . En , 

— g pp (m p c + — )x 



J p (t,m p c +j£-,z{). (29) 



Here, K n ~ 0.17 is the fraction of the proton kinetic energy 
that is transferred to the gamma rays or leptons. The num- 
ber of produced pions, fi, is a free parameter of the model 
that is fixed by requiring continuity of the gamma-ray emis- 
sivity at E-y = 100 GeV. 

The spectral intensity of gamma rays emitted from the 
jet can be obtained from 



7 7 (t, i? 7 



7r(zjtan£) n 



rhjZo dN~f(t, E n 



2 dNj(t,Ej,Zj) 



dE-, 



m p Vb 



dE-, 



(30) 
(31) 



We show the obtained result for gamma rays produced in 
the approaching jet in Fig. [4] 



3.4 High-energy neutrino emission from pp 
interactions 

Neutrinos are produced by the decay of the charged pions 
resulting from pp interactions, 



7T 
7T H 



e v u v e v u 



(32) 
(33) 



In the present context, the equipartition condition implies 
a strong magnetic field at the base of the jet. Hence, the 
charged mesons will undergo important losses due to syn- 
chrotron radiation. In order to find the local steady distri- 
bution of the parent charged pions at fixed E p , we have 
considered the following transport equation 



dN^(t,E-y,Zj) 
dK 



= 2 



s/E^ — m'ic 4 



dE^ 



(28) 



dKjE^F^Ey) JVCgQ 
dE n TTt(E n ) 



E^ ^ (E n )t pp . 



(34) 
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Here, t„ = t^ -^, with ri°' = 2.6x10 8 s, is the pion mean 
lifetime, b 7T (E 7T ) = (t^ c C^ir, Zj) + t^ iab (Zj)), and 



AaB^x" 
r(l 



+ 



1 

2x1 



1 



rx°>(l-x%) 



1 - 



%1Y Ep 



1/2 



(35) 



is the distribution of injected pions per pp collision 
|Kelner et all 120061 ), where a;^ = EV/Sp, B» = a + 0.25, 
a = 3.67 + 0.83L + 0.075L 2 , r = 2.6/ya, and a = 0.98/Va- 
The transport equation (1341) includes the effects of decays 
and energy loss of pions in the left member and the injec- 
tion of pions in the right side. 

The corresponding solution can be written as 



exp 



\K(E W )\ 

1 1 



+ ^log( *± 

b z E^ b z E> b'l & \E' 

a z / b z + a z E' 
of ° S \b x + aA 



dE', (36) 



where 



a z = - 



8n mecim-xC 2 ) 



2 v b r. 



The spectrum of high energy neutrinos from the direct 
decay of the steady distribution of pions is then 



F v (x,E v 



F K 



E„ 



,E P 



d.r 

X 



(37) 



where x = E v /E p and A = 0.427. 

As it is the case for the pions, the muons also undergo 
synchrotron and adiabatic losses. In this case, their mean 
lifetime is much longer (r^ = t^j^ with r^ ' = 2.2 x 10~ 6 
s). This implies that these leptons will lose most of their 
energy before decaying, especially at the inner parts of the 
jets, where most of the emission is originated. This can be 
seen in Fig. [5] where we show the decay rate t^ 1 and the 

, as a function of the muon 



loss rate iJ oss 



energy, for different values of Zj. According to this figure, 
muons with energies beyond 1 TeV will be present in the 
jets only at Zj > 10zo- Hence, since we are interested in 
neutrinos with energies E u > 1 TeV, the neutrino emission 
is attenuated due to the synchrotron losses in our model. 
The above equations show that to a good approximation 
we can safely neglect the contribution from muon decays at 
high energies. 

The neutrino emissivity can then be expressed as 



dN v {t,E v , 
dE v 



— U K — 



x 



F v ( x, — ) dx (38) 



(39) 



dN v {t,E v 
dE v 



The total neutrino spectral intensity emitted in the jet thus 




10° 10 1 10 2 10 3 10 4 10 5 10 6 10 7 10 8 
E u [GeV] 

Figure 5. Decay and loss rates for muons at different values of 




Figure 6. Spectral intensity of neutrinos produced in the ap- 
proaching jet as a function of the precessional phase and the neu- 
trino energy. 



reads 



rhjzo dN„(t,E v ) 
rripVb dE v 



(40) 



We show the result obtained for neutrinos produced in the 
approaching jet in Fig. [S] 



4 OPACITY TO GAMMA-RAY 
PROPAGATION 

The various mechanisms at work for gamma-ray absorp- 
ti on in the microquasa r SS433 have recently been studied 
in lRevnoso et all (|2008l ) . Absorption can occur via 77 inter- 
actions with ambient soft photons and by 7./V interactions 
with disk and star matter. 

In the first case, absorption occurs through interactions 
with low energy photons originated in the star and in the 
extended disk. These take place in the pair creation channel, 
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and the corresponding optical depth results from integration 
of 

dr 77 = (1 — e 7 • e p h)w p h(-E, f2')<7 77 dp 1 dE dcosd' d(f> (41) 

as described in lRevnoso et al] l|2008l ). Here, dp 1 is the differ- 
ential path followed by the gamma ray, E is the soft photon 
energy, e 7 is the unit vector in the direction of the gamma 
ray, and e p h = (sin 9' cos <f)' , sin 9' sin 4>' , cos 9') is the vector 
directed along the direction of the soft photons. The cross 
section for the process 77 — > e + e~ is 



u 77 (E^E) 



7rro 



(l-£ 2 )x 

- 2) + (3 - C 4 ) In 



i+i 



where t-q is the classical electron radius and 



£= 1- 



2(m e 



£ 7 £(1 



1/2 



(42) 



The radiation density of soft photons, in units 
cm^erg^sr- 1 , is n ph (E,Q') = 2E 2 {hc)^ (e E ' kT - l)" 1 
with T = 8500 K for the starlight photons and with 
T = 21000 K for the UV photons from the extended 
disk. The mid-IR emission is characterized by a radiation 
density n pb (E, fi') w F m d 2 / (hcE-Kr^ cos(0.62tt)), with 
Fir = 2.3 x lQ- 23 (A/Mm)-°- 6 for 2 /jm < A < 12 pxa 
l|Fuchs et al.ll2005l ). 

As for the absorption due to interactions with matter, 
the important effects are photopion production yN — > 7r'7 
and photopair production jN — > Ne + e~, where N repre- 
sents a nucleon. This last effect has not been previously 
considered, so it will be taken into account in the present 
work. We assume that the star has a matter density 



P*{r) = 



Ah 



ATvR^r' 



:Q(r - 



where r is the distance from the gamma-ray position to the 
center of the star. For the extended disk, we consider that 
matter density is given by 



,) = 



7j w Af2r 7 ' 



for 60° < 9z < 120°, where 9z is the polar angle in a coordi- 
nate system with its Z-axis directed along the approaching 
jet axis (for details, see Reynoso et al. 2008). The -yN con- 
tribution to the optical depth is 



r 7 jv(ij) = / <j 1 n 



_(p* +Pw) 



dp 1 



(43) 



where <t 7 jv = 0$ + a$ can be obtained from equations 



lp| and (fT6|). 

In Fig.[7]we show the total optical depth as a function of 
the precessional phase for different energies of gamma rays 
originated at the injection point 20 of the approaching jet. 
A very similar result is obtained the base of the receding 
jet, since zo is much smaller than any other size scale of 
the system. We clearly see the peaks of extreme absorption 
produced every time the star eclipses the emission region. 
The dependance of the total optical depth on the distance 
to the black hole Zj is also shown in Fig. [8] for gamma rays 
of energy E-, — 1 TeV coming from the approaching jet. 
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Figure 7. Total optical depth as a function of the precessional 
phase for gamma rays of different energies originated at the in- 
jection point zq of the approaching jet. 




Figure 8. Total optical depth as a function of the precessional 
phase and Zj for gamma-rays with Ej = 1 TeV originated in the 
approaching jet. 



5 NEUTRINO AND GAMMA-RAY 
DETECTABILITY 

The spectral intensities shown in Figs.|4]and|6]were obtained 
for g re i = 10 -4 without considering absorption. However, it 
is of course necessary to include the absorption effects to see 
how they affect the produced fluxes that may arrive to the 
Earth. 

The differential gamma-ray flux to be observed from 
each jet can be obtained as 

d$ 7 (/j,_E 7 ) _ 1 



dE-, 



47rd 2 



7r(zj tan£) n p x 



dN 1 {t,E 1 , zj) e - rtot (t,fi 7 , Zj ) 
dE-, 



dzi. 



(44) 



(45) 



We show in Fig. [9] the joint contribution of both jets 
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Figure 9. Differential gamma-ray flux arriving at Earth as a 
function of precessional phase and energy. 



Figure 10. Differential neutrino flux arriving at Earth as a func- 
tion of precessional phase and energy. 



to the differential gamma-ray flux considering g rc i = 10~ . 
Here, the absorption effects and the precessional phase be- 
havior have become manifest in the spectrum (c.f. Fig. 3}. 

As for neutrinos, although they undergo only weak in- 
teractions, we can estimate the corresponding neutrino op- 
tical depth using an expression analogous to (|43p . 



TV AT 



(3) 



{p* +Pw) , 
<7 i/JV dPy ■ 



(46) 



10 "'p 

Here the total u^N cross section can be approximated for 
E v > 1 TeV as a vN (E v ) » 10"" (B " } cm 2 , with 



au{E v ) = -38.42+1.46 log ( — ) -0.116 log 2 ' ' 



VGeV 



0.0041 log 



\GeV 
GeV 



The differential neutrino flux arriving to Earth can 
therefore be estimated as 



d<f> v {t,E v ) 
dE v 



87rd 2 



7r(zjtan£) n p x 



dN v (f E v ,z{) e _ TvN(ttEutZi) 
dE v 

Notice that an additional 1/2 factor has been put in order to 
take into account the reduction in the muon neutrino flux 
due to flavor oscillations over astrophysical distances (e.g. 
Athar et al. 2005). The result for the differential neutrino 
flux from the two jets with q Te \ = 10 -4 is shown in Fig. [10] 
as a function of precessional phase and neutrino energy. 

We can now integrate the differential fluxes in energy 
to appreciate the precessional dependance of the signals, 



It is interesting to explore the energy ranges that are ex- 
pected to be covered with different instruments. In the case 
of neutrinos with E v > 1 TeV, in view of possible detec- 
tion with IceCube, we obtain the flux shown in Fig. 1111 



where the contributions from the two jets appear sepa- 
rately. It is also included in that figure the value of the 
upper limit that can be extracted from AM ANDA-II data , 
$l lim) (£„ > 1 TeV) = 2.1 x 10 _11 cm _2 s _1 (|Harzenll2006h . 
which agrees with rece ntly published experimental results 
l|Achterberg et aHl2007r ) . The expected sensitivity for a km 3 
neutrino telescope such as IceCube, <I?! km3 '(i5„ > 1 TeV) w 
2 x 10~ 12 c m~ 2 s~ 1 for th r ee years of operation, i s also shown 
in Fi g. ITT1 (|Halzenll200"6l ; iDistefano et alj|2007l ; lAiello et al] 
l2007h . 

As for gamma rays, the sensitivity expected for GLAST 
is $^ LAST) ~ 6 x 10" 9 cm~ 2 s _1 at 100 MeV < £ 7 < 300 
GeVQ while for Cherenkov telescopes such as VERITAS 
and MAGIC II we consider a sensitivity ^( Cher ) ~ 9 x 
10~ 1 2 cm~ 2 s~ 1 for energies Ej > 100 GeV (Ca rmona et al.l 
2007). We show in Fig. [12] the integrated gamma-ray fluxes 
expected for q TC \ — 10 -4 in the mentioned energy ranges, as 
compared with the respective sensitivities. There, the con- 
tributions from the two jets are plotted separately. We also 
include in the lower panel the results corresponding to E-, > 
800 GeV, since they allow a comparison wit h the upper limit 
given by HEGRA for that energy range aronian et al.l 
l2005h . ^ lim) (S 7 > 0.8 TeV) = 8.93 x lO^cm-V 1 . 

According to our model, the most favorable range of 
precessional phases for the gamma-ray detection is found to 
be between ip < 0.09 and ip > 0.91, as it can be seen from 
Fig. 1121 This corresponds to the cases when the extended 
disk is so much open in our direction that the gamma rays 
emitted at the innermost regions can escape without having 
to travel through the thick extended disk, undergoing only 
77 absorption. 

However, since the observations by HEGRA were not 
perfor med exclusively at the d ates of favorable precessional 
phases l|Aharonian et al. 2005), we compare the value of the 
HEGRA cut with that of the averaged flux over one pre- 



3 See the official NASA GLAST webpage, 
http: / /glast. gsfc.nasa.gov/science/instruments/tablel-l.html 
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Figure 11. Neutrino fluxes arriving at Earth with E v > 1 TeV as 
a function of the precessional phase. The contributions of the two 
jets are shown separately, solid line for the approaching jet and 
dashed line for the receding one. The upper limit from AMANDA- 
II data and the expected km 3 sensitivity are shown in grey solid 
line and dotted line respectively. 



cession cycle, according to our model. Thus, the implied 
maximum value allowed by the HEGRA limit for the free 
parameter q Te \ is 



(max) 
ill 



10" 



<*7> 



2.9 x 10" 



which is below the maximum value allowed by the 
AMANDA II limit. 

Given that this parameter can be linearly factored out 
in our equations provided that q Ic i <C 1, which is guaran- 
teed for the allowed cases, we can plot the predictions for 
the different averaged fluxes as a function of this parameter 
and compare them with the expected sensitivities for the 
different instruments. 

This is done in Fig. 1131 where the mean neutrino flux as 
a function of q Ic i is shown in the upper panel as compared 
with the AMANDA-II upper limit and with the expected 
sensitivity for km 3 neutrino telescopes such as IceCube. The 
gamma-ray fluxes for 100 MeV < E n < 300 GeV and for 
E-f > 100 GeV averaged in the above mentioned range of 
favorable precessional phases are shown in the middle and 
lower panels of the same figure, as a function of q IC \. The ver- 
tical line indicates the bound q^^ derived from HEGRA 
observations, which is clearly more restrictive than the one 
implied by the AMANDA-II upper limit. 



6 DISCUSSION AND SUMMARY 

We have studied the high-energy emission originated in the 
dark jets of the microquasar SS433. A small fraction of its 
particle contents are relativistic protons that collide with 
the cold ions within the jets, producing gamma rays and 
neutrinos after pion decay processes. We found that up to 
distances ~ 10 12 cm from the black hole, protons with ener- 
gies below ~ 3 x 10 6 GeV will cool dominantly via pp interac- 
tions. The ratio of the power carried by relativistic protons 
to the total kinetic power of the jet, q rc i, was kept as a free 
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Figure 12. Gamma-ray fluxes arriving at Earth as a function of 
the precessional phase for 100 MeV < Ej < 300 GeV in the upper 
panel, for E-y > 100 in the middle panel, and for E 1 > 800 GeV 
in the lower panel. The contributions of the two jets are shown 
separately: solid line for the approaching jet and dashed line for 
the receding one. 



parameter of the model (for illustrative purposes we have 
used q rc i = 10~ 4 in the figures). 

We have calculated the spectrum of gamma rays and 
high-energy neutri nos based on the formulae given by 
iKelner eta l. (2006). We have considered the cooling of the 
charged pions and muons produced, and we have found that 
the high-energy neutrino emission is attenuated by syn- 
chrotron losses. Adding the contribution from both jets, 
we have obtained the total gamma-ray spectral intensity of 
SS433. The gamma radiation will be largely absorbed while 
leaving the inner regions the system by means of several 
processes. This will mainly occur through interaction with 
matter of the star and extended disk, leading to significant 
photopair production. UV and mid-IR photons originated 
in the extended disk are also expected to cause important 
absorption via 77 annihilations. 

The total optical depth is found to depend on the pre- 
cessional phase in such a way that when the approaching 
jet is pointing away from the Earth, at ip ~ 0.5, the ex- 
tended disk blocks the emitting region and the absorption is 
strongest. In particular, in the range of precessional phases 
between ip > 0.91 and tp < 0.09, the gamma rays originated 
at the base of the jets will travel to the Earth without having 
to pass through the equatorial disk. With the conservative 
assumption that this outflowing disk presents a large half 
opening angle a w = 30°, the mentioned range of favorable 
precessional phases gives a tot al of ~ 29 d for op timal de- 
tectability. Since, according to iGies et ail l|2002bl ). ip = 
occurred on 2002 June 5, it follows that the next upcoming 
opportunities to achieve detection will be centered around 
the following dates every 162 d: 2008 August 20, 2009 Jan- 
uary 29, 2009 July 10, etc. As mentioned above, the exact 
duration along which the favorable conditions may hold, de- 
pends on the half opening angle of the extended disk, which 
might be smaller than what was assumed. In that case, the 
observational window would be broader 
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Figure 13. Dependance on q Te \ of the different flux predictions: 
average neutrino flux {upper panel) and average gamma-ray flux 
for 100 MeV < Ej < 300 GeV (middle panel) and for E 1 > 100 
GeV (lower panel). The vertical line indicates the maximum pos- 
sible value of q rc i derived from HEGRA observations. The corre- 
sponding sensitivities are also show (see the text for details). 

The observations by HEGRA imply a maximum value 
for the free parameter q rc i in ~ 3 x 10~ 4 . Given the expected 
sensitivity of the next km 3 neutrino telescopes generation, 
it will be possible to test our model down to values ~ 5 x 
10~ 5 in three years of operation. An extended range of this 
parameter will be probed by gamma-ray observations with 
GLAST and Cherenkov telescopes, especially if performed 
on the favorable dates. 

We conclude that its dark jets can be the possible site 
for both gamma-ray and neutrino production in SS433. Since 
most of the high-energy flux is generated in the inner jets, 
gamma-ray absorption will make detection with Cherenkov 
telescopes like VERITAS and MAGIC-II difficult but not 
impossible if attempted at the favorable dates when the ap- 
proaching jet is closest to the line of sight. Actually, there 
are much better prospects for gamma-ray observation with 
GLAST and neutrino detection with IceCube also seems 
promising. The determination of the gamma-ray to neutrino 
flux ratio would allow to estimate q Te \ unambiguously, yield- 
ing crucial information about acceleration mechanisms tak- 
ing place in the jets. 
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